General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



X-640-7M23 

PREPRINT 


NASA I. X- 65^93 


CHARACTERISTICS AND CONSTANTS 
OF MOTION METHOD FOR 
COLLISIONAL KINETIC EQUATIONS 


FERDINAND F. CAP 


f 7. t **■ I •> ' 4 > l *01 

i\ 4 J • > ' 

(ACCESSI^pfUMBER) 

If AC tS) . \ 

T7r\X 


(THRU) 


(CODF.I 


^ (NASA CR OR TMX OR AD NUMBER) 


(CATEGORY) 


APRIL 1971 






0 6 


<• V 

V > 


GODDARD SPACE FLIGHT CENTER 

GREENBELT, MARYLAND 



CHARACTERISTICS ANT < CONSTANTS OF MOTION METHOD 


FOR COLLiISIONAL KINETIC jTQUATIONS 


Ferdinand F. Cap 

Theoretical Studies Branch 
Laboratory for Space Physics 
NASA/Goddard Space Flight Center 
Greenbelt, Maryland 20771* 


*On leave of absence from Innsbruck University, Innsbruck, Austria 


IWTROOUCTTON: 


The constants of motion n» thod* has been widely used to solve the 

+ 

Vlasov equation, and also for stability analysis . For this equation the 
nr- thod coincides with Lagrange's characteristic method: The Lagrange 

characteristic differential equations are identical with the equations of 
motion connected with the Vlasov equation. This is, however, not neces- 
sarily the case for other kinetic equations. 

In this paper we consider collisional kinetic equations and their 
characteristics. We also intend to use the constants of motion method to 
solve these equations. There are two main reasons for this investigation: 

1. This method of integrating kinetic equations is able to deliver 
exact nonlinear solutions (being of interest for nonlinear stability analysis, 
for nonlinear damping, etc.). 

2. Linearized solutions of the Vlasov equation breakdown when ^p 

3 o 

becomes the same order of magnitude as |p , Collisional effects may be 

3f 

proportional to fi or to and should therefore be included into a non- 

linear analysis (e.g. for a better matching together the linear and the non- 
linear solutions, mainly for big wave lengths and low frequencies) . 


THE COLLISIONAL KINETIC EQUATIONS 

The Lagrange characteristics of the Vlasov equation, coinciding with 
the respective equations of motions, read 


dx_+ du e a e 

at = u ' at = ‘ m e " se [u x g] 


a) 


In order to include collisional effects (encounters) due to long 

14 , 1 5 

range forces, we include a Iangevin term and obtain a new equation 


* References 1 through 6 
References 7 through 13 


of motion (Langevin equation) 


dx 

3t " 


u, 


du e + _ e 

5t m me 


[u x B] 



( 2 ) 


where v is an effective collision frequency. If such a dynamical fraction 
term is included, the Hamiltionian canonical equations of motions (which are 
used in the derivation of the Liouville equation resp the Vlasov equation) 
produce a kinetic equation 


!f+ <3»>f- £i + ^tSxSm u f-vf <3) 

of the Bhatnagar-Gross-Krook type lc . The Lagrange characteristics of (3) 
are, however, given by (1) and an additional equation 


1 df 

rat 


= v 


(4) 


We see that for the BGK equation (3) the Lagrange characteristic 
equations and the equations of motion are not the same. 

If one asks the question for which kinetic equation the characteristics 
are described by the equation of motion (2) and by (4) one gets another 
kinetic equation 


|| + (uV)f + (bV u )f = vf 


(5) 


An even more sophisticated collisional kinetic equation was given 

17 * 18/19 

bv Chandrasekhar and others 


3f 

3t 


+ (uvjf + 


(bV u )f = 3vf + qv^f 


( 6 ) 


where q is a constant. All these collisional kinetic equations (3) , (5) and 
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(6) nvay be derived also by various assumptions and neglections from the 
Fokker-Planck equation 2 .° 

The right-hand terms of the collisional kinetic equations are also 
of interest for particle belt physics: now the Liouville theorem docs no 

more forbid particle trapping since collisions within the field may change 
the particle properties 2 . 1 There is experimental evidence for such a change 
occurring in the Van Allen belt 21 and tnerefore in favor of the use of 
collisional equations. 


TRANSFORMATIONS OF THE KINETIC EQUATIONS 

Since we shall show later how to solve (6) , when we have a solution of 
(5), we first solve (5). We write 


f (x,u,t) = e g(x,u,t) 


and obtain from (5) 


|f + (uV)cr + (£v u )g = 0 


which is a Vlasov equation with a Langevin collision term. The Lagrange 
characteristic of (8) are given by the equations of motion (2) . So we 
may now solve the* collisional kinetic equations (5) , (6) and (8) and also 
(3) - by the constants of method. 

If ai...a 6 are six constants of motion of (2) 

a^x,u,t) = const (< 


so that 


u = u(ai . . .a 6 ,t) 
x = x(a 2 . . .a 6 ,t) 


( 10 ) 


% 
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is a solution of (2), then g(x,u,t) = g(ai...a6> is a solution of (8) 
since 


3a7 1 3t 1 + (uV)a i + ST ‘ 0 


Another useful transformation is the transformation to a wave frame 
moving with the constant velocity v. The variables in this wave frame may 


r = x -vt, w = u - v 


so that 


g (x,u,t) = 


.4 4 

F(r,w) 


presents a solution of (8) . The solution F is now determined by 


(WV r )F + (b\7 w )F = 0 


It is time independent if the fields E and B are time independent. 

2 o ”4 "4 # 

If one assumes F (r,w,t) (e.g. for special damping investigations) one 
obtains 


H + (WV r )F + (bVJF = 0 


CONSTANTS OF MOTION : 

In order to solve our collisional kinetic equation (8) we need constants 
of motion of the Langevin equations of motion. Multiplying (2) by u 


f 


we obtain 
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£ du _ e M dx _*dx 

u at ' s E(x > t) at ^ 


( 16 ) 


Since u = u (t) , x = x (t) , we may write u = u(x) and integrate 
2 

a (x,u,t) = 5 + | <J>(x,t) + vi^(x) (17) 

2 m 


where 


4 > = / E (x, t)dx, ip = /u(x)dx 


(18) 




The question if these integrals, exist depends on the special problem. 
The first equation of (2) gives (i = 1,2,3) 


b. (x,u,t) = t - /— 

1 A 


dx. 

l 


r. t - f: 


d XjL 
u. 


(19) 


2a-\^ -u^ -2$- .< p 


where i,k,l are cyclic permutations oi 1, 2, 3, x.^ = x, x 2 = y, etc. Now 
g (a,b^,b 2 ,b^) is a particular solution of ( 8 ) which can be proven directly 
by insertion. When treating special problems, two further constants of 
motion may be found. 

TOE CHANDRASEKHAR EQUATION : 

In treating ( 6 ) we follow ideas which were discussed by Chandrasekhar ! 7 ' 1 8 ' 1 9 
For simplification we discuss only the one-dimensional case 

3f . „ 3f . h 3f _ vr _ 9 2f ,„ m 

at 3x b au ^ gu 7 

Using (7) we obtain 


!f +u £ + b £=*!£ 


( 21 ) 




Introducing g(a,b,t) into (21) where a and b from (17> to (22) respect- 


ively we get 



Putting u 2 (t) = G(t) we obtain the particular solution 
g(a,t) = exp[q/G (t)dt-qt-a) 

The methods proposed in this note shall be applied to forthcoming 


(23) 


(24) 


papers to collisional nonlinear electrostatic modes, to collision.?! non- 
linear Landau damping, to collisional nonlinear damping of large amplitude 
whistler waves, etc. 
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